FOREWORD 



The first attempt to apply a mathematical framework to thermocapillary flows is attributed to the article written 
in 1959 by N.O. Young, J.S. Goldstein and M.J. Block [YGB]. This classic work brought interface driven flows 
into attention of the Eastern hydrodynamic community particularly focusing on the surface flows induced by the 
temperature inhomogeneities. However, there is another work, which is currently hidden from worldwide scientific 
attention due to the fact that it was written in Russian many years ago. In 1956, in Zhurnal Fizicheskoi Khimii (a 
journal of the former USSR, nowadays it is known as Russian Journal of Physical Chemistry A) A.I. Fedosov published 
an article entitled 'Thermocapillary motion', where he introduced a mathematical description of the thermocapillary 
effect, sequentially considering two problems: the motion of a flat liquid layer and the motion of a spherical non- 
deformed drop without gravity. After thorough historical research, and with much help from different people we 
found that this result was obtained by Fedosov before the year 1948, in his doctorate thesis (under the supervision 
of Benjamin Lcvich). However, the story is even more curious: in 1944, Lev Landau and Evgcny Lifshitz published 
[LL] the most general form of the boundary conditions for the liquid-liquid interfaces which ultimately leads to any 
cause of surface-driven motion, including the thermocapillary one. Below, we present an English translation of the 
Fedosov article, temporarily leaving aside the chronology of the scientific inputs of Fedosov, Levich and Landau, and 
reserving to return to this story in the near future. 

In a course of this translation, a substantial amount of historical work was done. We deeply appreciate Sergey 
A. Fedosov, Alexey V. Belyaev, and Igor Yu. Makarikhin for their generous help in obtaining the important docu- 
ments. We are also indebted to Michael Kopf for critical reading of this translation. 

[YGB] N.O. Young, J.S. Goldstein, and M. J. Block, The Motion of Bubbles in a Vertical Temperature Gradient, 
J. Fluid Mech. 6, 350-356 (1959). 

[LL] L.D. Landau and E.M. Lifshitz, Mechanics of Continuous Media (OGIZ, Moscow, 1944) (in Russian). 
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Mechanical motion (convection) appears in a nonuniformly heated liquid because its density de- 
pends on temperature 1 '. This flow causes the liquid to self-mix unifying the temperature field. The 
convective flows are addressed in a series of papers recently published [1]. When the liquid has an 
interface, besides the convection, a different kind of a motion can appear. The origin of this motion 
is a gradient of the surface tension 2 '. Similarly to the electrocapillary motion, one can call this new 
motion a thermocapillary motion. In the present paper we consider two classes of thermocapillary 
motion: motion of a liquid in an open container and motion of a liquid drop suspended in another 
liquid. 



Motion of liquid in a flat open container when a 
temperature gradient is applied along the liquid 
interface 



Considering the motion of liquid in a flat container let 
us assume that width and length of the container are suf- 
ficiently larger than its depth. This assumption reduces 
considerably the complexity of equations of liquid mo- 
tion. The component momentum equations for steady 
motion of a viscous incompressible liquid written in the 
Cartesian coordinates are 
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Here F is an external body force 3 ) . 

Let us introduce an origin of coordinate system on the 
interface of liquid. The x- and y-axes are applied along 
the container length and width, respectively, the z-axis is 
directed from the liquid interface down to the bottom of 
the container. We assume a constant temperature gradi- 
ent is applied along the x-axis, giving us v y = 0. 

We will also neglect the convective flow because the 
depth of the container is smaller than its length. Fea- 
sibility of this assumption and the limitations which it 
causes for the theory will be addressed below. Now, 



1 ) This is true with gravity applied (KM & SB). 

2 ' Currently, this kind of motion is known as Marangoni convection 

in contrast to gravitational Raylcigh-Benard convection (KM & 

SB). 

3 ) v = fJ,fp, where u, p,, and p are kinematic viscosity, dynamic 
viscosity, and density of the liquid, respectively (KM & SB). 



since we disregarded the convective motion, then far away 
from the container walls v z = 0. Further, the term 
dv x /dx ~ v x /l (I is the container length) is much smaller 
than dv x /dz ~ v x /h {h is a depth of the container) and 
can be omitted, too. As a result of this simplifications, 
the system (1) is reduced to only two equations. The first 
one is 
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Liquid in the container is dragged by the moving surface 
layer, then it returns backward due to the container walls 
causing the appearance of a pressure drop along the x- 
axis. The second equation is 
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After integrating this equation the solution is 

Pi{x) 
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pgz. 
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(4) 



Now we will define the boundary conditions for Eq. (2). 
First, on a solid boundary the liquid velocity vanishes, it 
gives us a first condition 

v x (z = h) = 0. (5) 

Second, on a free liquid interface the components of a 
viscous-stress tensor are continuous. Since v y = v z = 0, 
then only the one tensor component p xz is essential for 
us while the others are either equal to zero or are not 
important at all. The continuity of p xz results to the 
second boundary condition 



da_ 

dx 



dv x 
dz 



0. 



(6) 



where a is the surface tension. Besides the boundary con- 
ditions, the solution of Eq. (2) should meet an additional 



condition representing the fact that the average velocity 
over the container cross section equals to zero: 



1 

h Jo 



v x dz = . 



(7) 



With due account of relation (4), equation (2) is easily 
integrated. Indeed, substituting Eq. (4) in Eq. (2) we 
have 

dpi(x) _ d 2 v x 



dx 



dz 2 



Since P\(x) does not depend on the coordinate z, the 
solution of Eq. (2) is 



v x = A + Ba 



1 dpi(x) 
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Applying the boundary conditions we have 
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and finally the expression for the velocity is 

1 da 1 d Pl (x) 2 2 

v x = T1 -(h -z)- ih -z). 

pi ax Ip, dx 

Applying this formula to condition (7) gives 

dp\(x) 3 da 
dx 2h dx ' 
and for the pressure in the liquid we have 

P=Po + -^W(x)-a(0)}+pgz, (8) 

where the constant po cannot be found because the pres- 
sure in liquid is always defined to the extent of an additive 
constant. Finally, it gives us the velocity function in the 
liquid 
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The profile of the velocity field is shown in Fig. 1. 

One can see from Eq. (10) that velocity of the liquid 
reaches its maximum on the liquid interface 



h da dT 
Ap dT dx 



(11) 



Since da/dT < 0, therefore along the interface, the ve- 
locity direction is opposite to the temperature gradient. 




FIG. 1. The profile of the velocity field for thermocapillary 
motion in a flat long container. 



If one compares the expression for the maximal ve- 
locity of the thermocapillary flow 4 ) against the expres- 
sion for the characteristic velocity of the convective flow, 
introduced below, then one may think that the convec- 
tive flow could be simply disregarded for the condition of 
large container depth h. However, it is not the case. 
In fact, formula (11) cannot be used for any size of 
the container. Indeed, in the momentum equations we 
have neglected the term v x (dv x / dx) ~ v 2 /l comparing 
to u(d 2 v x / dz 2 ) ~ v\v x \/h 2 . In other words, we assumed 
that v 2 /l <C u\v x \/K 2 or, equally, h 2 <^ vlj\v x \. Substi- 
tuting Eq. (11) into the last expression, we have a validity 
criterion for the developed theory: 
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When this inequality is not valid, our theory is not valid 
too. Thus, we cannot argue that if the depth, h. of the 
container is large enough, then velocity of the thermo- 
capillary flow overcomes velocity of the convective flow. 
One would also think that with unlimited growth of both 
sizes h and I, the thermocapillary velocity could exceed 
its convective value, however it is not the case as well. 
There are two restrictions regarding this way. 

a) The temperature gradient is constant for any size 
of the container. In this case, the temperature difference 
AT over the container is AT = T - T = VTl - /i 3 , 
and the convective velocity 5 - 1 grows much faster than the 
thermocapillary one because the former depends on tem- 
perature difference AT and not on temperature gradient. 

b) The temperature difference is constant. The ex- 
pression for the maximal thermocapillary velocity can be 
rewritten as 
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Consider the estimate most favorable for the thermocap- 



4 > see Eq. (11) (KM & SB). 

5 > see Eq. (*) below (KM & SB). 
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illary flow 6 ' 



Then wc have 
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Thus, the maximal velocity increases with growth of 
the temperature difference and temperature coefficient of 
surface tension and decreases with the growth of viscosity 
and density. This is obvious. The most important result 
is the dependence of maximal velocity on the container 
size I. It is seen that under fixed temperature difference 
the maximal velocity of thcrmocapillary flow decreases 
with growth of /. 

Assuming the following values for the parameters /i = 
0.01 P, dT/dx = 0.1 K/cm, h = 3 cm, dajdT = 
—0.15 erg/(cm 2 K), involved in the thermocapillary ve- 
locity expression, we calculate the maximal velocity value 
v x ,max = 1-1 cm/s 7) . 



Let us compare this estimate with the value of the con- 
vective velocity at the same conditions. The expression 
for the convective velocity has a form [1] 8 ) 



where 6 is the boundary layer thickness, x < S and 



(*) 



5.17i/ 



5 = 3.93 (1" 1 



20 
21 



-1/2 



9(Ti - T ) 



1/2 



,1/2 



20 
21 



1/4 r 



9(Ti - T Q 



-1/4 



,1/4 



6 ' a limit of validity for the theory (KM & SB). 
7 ' The aforementioned part of the article is identical to the text in 
the PhD thesis [3] of Fedosov and to the part of the chapter 'Ther- 
mocapillary motion' of the book [4]. Fedosov's next consideration 
(the text in between two markers T and A) of the convective contri- 
bution is not applicable to the given problem. As it was shown by 
Birikh [5] , for a thin layer the convective contribution can be found 
analytically in the same way as the thermocapillary one. The later 
prevails for layer with thickness < 1 ~2 mm, sec details in [5] (KM 
& SB). 

s ' The problem considered in Rcf. [1] and cited here is the problem 
of a convective flow near the hot vertical wall. The geometry of 
the vertical wall problem is essentially different from the geome- 
try of the horizontal layer studied here. This mismatch between 
the geometries of the vertical walls and horizontal layer results to 
the incorrect estimate of the convective flow J CO nv (see the next 
footnote) (KM & SB). 



Here g is the gravity, T\ and To are the temperatures of 
the wall and far from the wall, respectively, v is the kine- 
matic viscosity, 7 is the thermal diffusivity, z is the non- 
trivial coordinate of a flat boundary layer. Since Ref. [1] 
does not provide an explicit expression for the convective 
velocity outside the boundary layer, we have to compare 
liquid fluxes rather than velocities: 

,5 x 
Jcowv pi i^QQYivdx — —u z 8p . 
Jo 6 

Substituting here u z and 5 we find 9 ) 
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The corresponding relation for the thermocapillary flux 
equals 



Jth- 
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Using the values of parameters mentioned above and 
p = 1 g/cm 3 , Ti - T = 1 K, T = 300 K, we obtain the 
following estimates for both flows: J C onv = 0.23 g/(cm-s) 
and Jth-cap = 0.5 g/(cm-s). Thus, the thermocapillary 
mass flux proves to be prevailing over the convective one. 
Note, the given estimates are true for fairly large temper- 
ature gradients 10 ^. If the temperature gradient was of the 
order of 0.01 K/cm n \ our conclusion would be changed 
to the opposite one and now the convective flow would 
dominate. Thus, one can conclude that the thermocap- 
illary motion is negligible in all practically interesting 
situations. However, this is not true for some experi- 
ments where the heating of liquid is caused by applying 
the light on its surface. In this case a liquid is heated in 
a very thin layer near the free surface. For the consid- 
ered termocapillary problem, this 'skin layer' doest not 
mean anything because the thermocapillary velocity de- 
pends on the total thickness of the liquid layer. For the 
convective flow, however, this 'skin layer' has a principal 
meaning because it is a variable in the velocity expres- 
sion. For the lower values of the 'skin layer' ~ -L of the 
total thickness of the liquid layer, the thcrmocapillary 
flow dominates again. At lower values of the 'skin layer' 
just the convective flow becomes negligible. 



9 ) The relevant consideration of the convective term has been done 
in [5]. There was found that the convective flux (or velocity) 
J co „ ~ Gr ~ (Gr is Grashof number), whereas the thcrmo- 
capillary flux Jth-cap ~ Mg ~ h 2 (Mg is Marangoni number). In 
the case of sufficiently thin layer, the thcrmocapillary flow always 
prevails (see details in [5]). Fedosov's estimate proves to be a fourth 
root of the correct one, J conv ~ Gr 1 / 4 (KM & SB). 

10 ) dT/dx = 0.1 K/cm (KM & SB). 

n ) e.g., more realistic for the long and shallow containers than the 
one above (KM & SB). 
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Motion of a drop in viscous medium due to 
temperature gradient 

Let us consider now the problem of motion in a temper- 
ature field of a liquid drop suspended in another liquid. 
Because of the temperature difference at different points 
of the drop surface, the later cannot remain at rest. In- 
stead, it will move from the warmer regions with the 
lower surface tension toward the colder regions with the 
higher surface tension. In the same moment the moving 
drop will drag the surrounding liquid medium by apply- 
ing some force to it. There will also be an equal force 
in opposite direction which the surrounding medium will 
apply onto the drop. This reactive force will cause a drop 
movement in the direction of the temperature gradient. 

Let us estimate the order of the thermocapillary ve- 
locity of the drop. The characteristic force applied per 
unit length of the nonuniformly heated surface is of the 
order of |^VT|. This for ce is balanced by the viscous 
forces 12 ^: 



a a 
Therefore we obtain 
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The exact value of the thermocapillary velocity can be 
found by solving the hydrodynamic equations with the 
appropriate boundary conditions. We will use the hy- 
drodynamic equations within the Stokes approximation, 
because in all practical cases both the thermocapillary 
velocity and the Reynolds number are small 

Vp = /iAv , Vp' = fi'Av' , 
divv = 0, divv' = 0. (13) 

Here and below the variables with primes denote the 
liquid domain of the drop, while the variables without 
primes relate to the surrounding liquid. Let us consider 
the problem in the polar coordinate system attached to 
the moving drop and with the origin in the drop center. 
We chose the polar axes along the temperature gradient. 
The velocity is defined to be positive when it is directed 
along the polar axis. The boundary conditions far from 
the drop are 

(14) 
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At the drop surface the following components of the stress 
tensor are continuous 13 ) 

(.Prr Prr)\r—a Prr\r— & 7 



12 ' fi and \J , u and u' are the dynamic viscosities and liquid veloc- 
ities inside a drop and in the surrounding liquid, respectively, a is 
a drop radius (KM & SB). 

13 ) The bulk and surface components of the stress tensor are shown 
here by upper indexes '0' and '<r', respectively (KM & SB). 
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In addition, at the drop surface the following conditions 
are true 



0. 
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Here v r , vg, v' r , v' g are the corresponding components of 
the velocity field. The stress tensor components are 
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Expressions for the variables with primes, p$^, are 
analogous to those for p® r and p® s with the obvious sub- 
stitutions of unprimcd functions by their primed counter- 
parts. Because of the uniaxial symmetry of our problem, 
the velocity component v v = and all the variables of 
the problem arc independent on the azimuthal angle ip. 
The surface tension is a part of boundary conditions and 
depends on the distribution of temperature on the drop 
interface. 

We assume that the motion of the drop does not af- 
fect the distribution of temperature in the surrounding 
liquid, because the drop moves very slow. The temper- 
ature inside the drop is as if the drop had stayed at its 
particular location long enough. In addition, we assume 
the temperature gradient to be constant. Therefore the 
temperature distribution in both the drop and the sur- 
rounding liquid is described by the expression 



T = T r 



|VT|rcos6>. 



(18) 



where, T c is a temperature in the center of the drop. T c 
is a function of time 



T c = T + \VT\ut . 



(19) 



Here To is temperature in the center of the drop for initial 
drop position and t is time. We also assume that the 
motion does not affect the drop spherical shape. Later, 
we will consider the feasibility of this assumption as well 
as its limits. 

When the above assumptions are true, one can take 
the first spherical harmonic functions as a solution for 
the hydrodynamic equations [2]. 

The gradient of the surface tension caused by the in- 
homogeneous temperature on the interface is an origin of 
the motion on the drop surface. One can find the surface 
tension as a function of the angle 9 using the following 
expression 
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As a first approximation, one can assume that the deriva- 
tive da/dT does not depend on temperature and there- 
fore on angle 8. In this case, taking into account Eq. (18) 
one finds the surface tension 
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Now we insert the last expression to the boundary con- 
ditions (16) getting seven equations to obtain the eight 
constants ao, bo, ai, b\, a 2 , b 2 , &3, u 14 ' ) ■ As usual, one 
constant cannot be defined, so let us set ao = po, where 
po is a pressure in the center of the drop. The other 
constants are: 
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The motion of liquid outside and inside the drop are de- 
scribed, correspondingly, by the following equations: 
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Velocity u (this is velocity of the liquid at infinity) is op- 
posite to the temperature gradient because of the rela- 
tion da/dT < . Velocity of the drop is — u, it is aligned 
along the temperature gradient. Thus, the drop moves 
from the cooler layers of liquid to the warmer ones. It 
is apparent that this motion promotes the temperature 
leveling-off. 

Let us compare the thermocapillary velocity of the 
drop with a maximal velocity of the thermocapillary flow 
in the container 
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Now, if one uses the characteristic values for the pa- 
rameters fj, = 0.01 P, n' = 0.5 P, a = 0.3 mm, |VT| = 



0. 1 K/cm, da/dT = -0.015 erg/(cm 2 -K) in the drop ve- 
locity expression, then one can estimate the drop velocity 
u = 2- 10~ 4 cm/s 15 '. 

We have mentioned above that assumption of the con- 
stant spherical shape of the drop should constrain our 
theory. Let us consider this constraint now. We believe 
that the drop keeps its shape nearly spherical when the 
surface tension on the drop interface is nearly constant, 

1. e., \a — a^/ 2 \ 0V/2j where ov/2 is the surface tension 
on the drop's equator and a is the surface tension in the 
any other point of the drop interface. Assuming this wc 
can write 

(22) 



(da/dT)\VT\a/a^ /2 < 1 



To make clear in which parameter regions our theory is 
applicable, let us consider a numerical example. The 
equator's surface tension depends on the equator's tem- 
perature, however, it is not significant regarding our es- 
timate of the theory's validity. Therefore for the values 
of the parameters given above we have |VT|a <C 50 K . 
This inequality requires the temperature change along 
the drop perimeter to be less than 50 K. It is clear 
that this inequality is almost always true. Indeed, for 
|VT| = 0.1 K/cm we should require the following in- 
equality a 5 m to be true. It is obvious that all drops 
in real experiments satisfy this condition. 

One can show that the thermocapillary drop velocity 
can overcome (and sometimes, for some cases does so 
significantly) the electrocapillary drop velocity for weakly 
conductive drops. 

I would like to thank V.G. Levich for giving this prob- 
lem to my consideration and for his constant advice dur- 
ing my work. 

CONCLUSIONS 

1 . Two cases of the flow of liquid caused by the tem- 
perature gradient in the interfacial layer (thermocapillary 
flow) are considered. 2. Thermocapillary flow can over- 
come convective flow when the interfacial heating layer 
is reasonably small. 3. The limits of the theory's appli- 
cability arc discussed. 
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